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We investigate the excitation spectrum of a two-dimensional resonating valence bond (RVB) 
state. Treating the 7r-flux phase with antiferromagnetic correlations as a variational ground state, 
we recover the long wavelength magnon as an "RVB exciton". However, we find that this excitation 
does not exhaust the entire spectral weight and the high energy spectrum is dominated by fermionic 
excitations. The latter can be observed directly by inelastic neutron scattering and we predict 
their characteristic energy scales along different high symmetry directions in the magnetic Brillouin 
zone. We also interpret experimental results on two magnon Raman scattering and mid-infrared 
absorption within this scenario. 
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The undoped high T c cuprates, such as La-iCuO^ etc., 
|J and compounds like Cu(DC0 2 )2 -^D 2 (CFTD) §|] 
are ideal realizations of two-dimensional (2D) spin 1/2 
Heisenberg antiferromagnets on the square-lattice com- 
posed of Cu ions and, consequently, are interesting sys- 
tems to study. The antiferromagnetism that is observed 
in the undoped cuprates plays a crucial role in some the- 
ories of the metallic state of doped Mott insulators. For 
these reasons, it is important to understand the nature 
of magnetic excitations in these materials. A direct way 
to probe the spectrum of magnetic excitations is inelas- 
tic neutron scattering (INS). The presence of long range 
antiferromagnetic order below the Neel temperature was 
established in the insulating cuprates by neutron scat- 
tering measurements 0. Well defined peaks, identified 
as magnons, have been observed to disperse along the 
magnetic zone boundary in different compounds 
Since then, a spin wave dispersion throughout the en- 
tire Brillouin zone in La.2CuOi, has been extracted from 
INS measurements j^] . The antiferromagnetic correlation 
length diverges exponentially as the temperature is low- 
ered, indicative of long range order at T = @||. This 
is in agreement with the analysis of the 2D antiferromag- 
netic Heisenberg model, with superexchange coupling J 
between adjacent spins, in terms of the quantum non- 
linear (7 model |(|, and other approaches based on spin 
waves |Q. These results, concerning essentially the low 
energy regime, are construed as evidence favoring the hy- 
pothesis that the excitation spectrum in these materials 
can be described adequately by renormalized spin wave 
theory f§. 

On the other hand, optical probes such as Raman and 
mid-infrared absorption spectroscopy highlight the inad- 
equacy of the spin wave hypothesis. Typically, the po- 
sitions of primary peaks observed in these experiments 
coincide reasonably well with calculations based on spin 
wave theory, but the observed line shapes cannot be ex- 



plained, owing to the presence of finite spectral weight 
at high energies. For instance, the B\ g shift in Raman 
scattering experiments spans a rather broad range of en- 
ergies in insulating cuprates Q. Recent measurements 
on insulating ST1CUO2CI2 as well as in YBaiCuzO§.\ 
show very broad line shapes, spectral weight at higher 
energies, and a broad feature around 4 J p0| . Optical 
absorption spectra show similar features. Experiments 
on insulating cuprates [ pdp^ j show a primary peak iden- 
tified with bimagnon plus phonon absorption ]l3| . But 
as in Raman scattering, the spectrum has a long tail ex- 
tending to 8000 cm -1 in La^CuO^ and 6000 cmT 1 in 
Y Ba2CuzO§. By analyzing the lineshape and the high 
energy spectra, Griininger and coworkers have presented 
strong experimental evidence against the spin wave hy- 
pothesis for the insulating cuprates [ p2[ . 

The two classes of experiments described in the preced- 
ing paragraphs do not necessarily contradict each other. 
Though the dispersing peaks observed in INS measure- 
ments can be interpreted as spin waves, an analysis of 
the spectral weight shows that long range order and spin 
waves can account only for about 50% of the observed 
spectrum |lj] . These results indicate the presence of ex- 
citations beyond the one magnon mode. Thus, taken in 
conjunction, INS and optical spectroscopy suggest two 
possible theoretical approaches. One is to start from 
spin wave theory, and look for consistent explanations 
for the experimental results outlined above. This is not 
an easy task, as spin wave theory is an effective theory 
of long wavelength excitations. The other approach is to 
postulate a new set of excitations that offers a natural 
explanation of the high energy spectra, and describing 
the magnons in terms of this new basis. It is the latter 
approach that we shall pursue. 

In this paper, we consider the flux-phase resonating 
valence bond (RVB) theory, treated within the random 
phase approximation (RPA). Gapless spin wave excita- 
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tions are recovered as a Goldstone mode. Besides this, 
the RVB theory also makes specific predictions about 
novel excitations- spin 1/2 flux fermions or spinons, and 
their direct observation. We predict the characteristic en- 
ergy scales of these excitations along various high symme- 
try directions and show how their presence can naturally 
account for the spectral weight seen at high energies in 
optical probes. 

Our starting point is the 7r-flux state with staggered 
magnetization as a variational parameter fl5|| . For 
the case of the 2D Heisenberg antiferromagnet, the 
above state yields the best variational energy when the 
Gutzwiller constraint is enforced exactly . We assume 
that the excited states of the flux fermion spectrum are 
also good variational states for the spin excitations of 2D 
Mott insulators. The 7r-flux state with staggered magne- 
tization m is obtained in the mean-field approximation 
of the Hamiltonian 

' H = -' I f E e^iflf^ + h.c^ + V^n^ , 

(1) 

which is an extension of the RVB mean field theory 
fj"7f . The lattice fermions pick up a phase $n=7r 
on hopping around an elementary plaquettc, and also 
experience an on-site potential V, which is included 
to induce a spin density wave (SDW) coexisting with 
the flux order. The excitation spectrum is given by 
Ek=J c ff yjcos 2 k x + cos 2 k y + (to/2) 2 where to is deter- 
mined self-consistently as V~ 1 =N~ 1 ^2' E^ 1 , N being 
the number of lattice sites and the summation being over 
momenta in the magnetic Brillouin zone (MBZ). 

The above Hamiltonian, with J g = J, was first consid- 
ered by Hsu p5| . Invoking a Gutzwiller approximation, 
Hsu determined the projected variational energy and the 
sublattice magnetization as a function of the SDW mass 
parameter to, induced by the potential V. The analysis of 
the excitation spectrum, however, is complicated by the 
Gutzwiller approximation. Owing to this difficulty, Hsu 
could only obtain the poles of the particle-hole Green's 
function (S z = ±1 excitations) and not the complete dy- 
namical spin susceptibility. To obtain the latter, we find 
it convenient to follow an equivalent approach, proposed 
first by Laughlin In this scheme, a spinon pair is 

created as a projected particle-hole excitation and the 
spinons constituting the pair experience an on-site repul- 
sion, favoring the formation of an SDW. This interaction 
also enhances the exchange integral J to J c g. The quan- 
tities J e ff and to are determined self consistently to be 
1.5J and 0.5, respectively [Q. The Gutzwiller approxi- 
mation technique employed by Hsu gives similar results. 

We now obtain the dynamical spin susceptibility from 
Eq.(|l]) using RPA. The calculation parallels that of Schri- 
effer et at, on the SDW instability of the 2D Hubbard 
model p9[ . We obtain, for the (transverse) spin suscep- 
tibility, 5(q, uj) = — Imx + ~(q, oj), where 



and 

x +-(q, w ) = --L^' 

k 

^ cos k x cos{k x + q x ) + cos k y cos(k y + q y ) — (to/2) 2 

-E-k-E-k+q 



l uj-(E k + E k+cl )+iS uj + (£ k + £ k+q ) -iS 1 

Here, \q~ is defined as the time-ordered product 
{TS + S~) with respect to the SDW ground state. 
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FIG. 1. Spin susceptibility S^q,^) as a function of energy 
uj for various q's in the direction of (a) to (7r,0), (b) 

to (n /2,n /2), and (c) (^,0) to (n /2,n /2) along the zone 
boundary. The results are obtained using Eq.(|2|) with 5=0.001 
and 3,000 q points in the MBZ. The energy is in units of 1.5J. 
The small vertical arrows mark the positions of the magnon 
peaks. 

At the magnetic wave vector q=Q=(7r,7r), x(q, oj) is 
singular when u>=0. The spectrum is gapless, which re- 
sult is guaranteed by the self-consistency of the RPA. 
Since the gapless mode is a bound state between particle- 
hole pairs, it is interpreted as an "RVB exciton" . The 
dispersion of this mode can be obtained analytically. 
Evaluating lim q ^Q for small oj, we obtain w(q) « 
1.14J e //|q|. For arbitrary q, Eq.(g) has to be evalu- 
ated numerically. The results for 5*(q, u>) are shown in 
Fig. |l|. For each q, we see a strong peak at low en- 
ergy. This magnon (or RVB exciton) peak disperses as q 
varies. By tracking the position of the intense peaks at 
low energy, we obtain the magnon dispersion for the en- 
tire MBZ (see Fig. |2j) , recovering Hsu's results [ fL5| . Near 
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the zone boundary, the dispersion deviates from that of 
linear spin waves. 

It is important to realize that our approach is much 
more than just a complicated way of describing magnons. 
This is borne out by our explicit evaluation of Sfq, u>). 
For example, consider S(q, uj), at (tt,w) in Fig. |[ The 
intensity drops at higher energy, and grows again at 
lu ~ 0.5J o ff. The intensity at higher energies is inter- 
preted as the manifestation of the particle-hole (spinon) 
excitations. Within this interpretation, we can estimate 
the characteristic energy at which the continuum be- 
comes visible. Consider a spin-flip excitation with mo- 
mentum q as a convolution of particle- hole excitations. 
Then, its energy w(q) is given by w(q) = -Eq-k + E^, 
where we have ignored the residual interaction. For 

q=Q, uj(Q)=2J c g^Jcosk^, + cosfc 2 + (to/2) 2 . Its mini- 
mum, occurring when the two fermions are created at the 
Dirac point (tt/2, n/2), is exactly mJ c s, as seen in Fig. []]. 
This analysis shows that the continuum emerging for en- 
ergies 0.5 J c ff and above can be interpreted as arising from 
a convolution of free flux fermions or spinons. The line 
shape is, however, determined by the interaction. Due 
to the phase space available for convolution, the onset of 
the continuum is a strong function of momentum. Hence, 
as seen in Fig. [I], the magnon peak as well as the lower 
boundary of the flux-fermion continuum shift in energy, 
as q is varied. As q — > (n, 0), the onset of the contin- 
uum shifts towards lower energies. At (71", 0), there is no 
easy way to isolate the magnon and spinon contributions 
pofl . In Fig. 0(b), we display the behavior of 5(q, ui) 
along the direction (tt,tt) — > (tt/2,tt/2). In this case, the 
lower boundary of the fermion continuum shifts towards 
higher energies. At (ir/2,ir/2), we find the onset of the 
continuum around 1.7J e ff. As before, this value can be 
understood by considering the spin-flip excitation with 
momentum (7r/2,7r/2) as a particle-hole pair. A similar 
analysis as in the previous case yields a minimum energy, 
J cff (m + 2- v /2 + (to/2) 2 )/2 at the Dirac point (vr/2,7r/2). 
For m = 0.5, this turns out to be 1.7J c ff, confirming our 
interpretation. 

From our results for the magnon peak and the fermion 
continuum, we obtain the spectrum of excitations over 
the whole MBZ, which is shown in Fig. [| Note that by 
putting = in Eq.(|l|), results from the weak coupling 
approach to the Heisenberg antiferromagnet jl9| can be 
obtained. We have verified that the excitation spectrum, 
in this case, is distinct from ours P0| . Thus, our calcula- 
tions are specific to the existence of a 7r-flux phase p0| . 
Our results for small q ~ Q are consistent with calcula- 
tions based on the bosonic RVB scheme |2^|. At higher 
energies, it has been shown that the bosonic RVB 
state can support vortex like topological excitations that 
are fermionic. This agrees with our picture of fermionic 
excitations at high energies. The shaded region in Fig. ^ 
depicts these excitations and can, in principle, be ob- 
served directly by INS. At the magnetic wave vector Q, 
there is a well defined gap between the strong magnon 



peak and the spinon continuum, which may facilitate de- 
tection of the continuum along this direction. For typ- 
ical values of the superexchange coupling J in undoped 
cuprates, this gap is of the order of 100 meV. Recently, 
Coldea et al. B have reported results from high-energy 
(0.1-0.5 eV) neutron scattering in La^CuO^. Given the 
possibility of studying high energy spin excitations with 
enhanced resolution, INS may be a good probe for the 
direct observation of spinon excitations, should they ex- 
ist. It is therefore extremely interesting to examine, both 
theoretically and experimentally, whether our predictions 
are realized in this measurement. Experiments show that 
as q varies from (jr, 0) to (tt/2,tt/2), the spectral weight 
at high energies decreases. This is certainly consistent 
with our calculation. A detailed comparison incorporat- 
ing appropriate structure factors for neutron scattering 
is forthcoming 20 1. 
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FIG. 2. Spectrum along high symmetry directions in the 
MBZ deduced from spin susceptibility intensities S(q,lj). 
Each (gray) diamond marks the 5(q,u;) which has in- 
tensity larger than 1.5 (arbitrary unit), as chosen here. 
The low-energy branch represents the magnon dispersion. 
The broad shaded region corresponds to the spinon contin- 
uum. The solid line represent s the spin-wave dispersion, 
2Z C J^J 1 — [(cos k x + cos k y )/2] 2 with Z c ~ 1.18 the renor- 
malization factor. 

Let us now turn our attention to the two magnon ex- 
periments. As mentioned earlier, these experiments are 
characterized by a broad spectral distribution of inten- 
sity and secondary peaks. Both these features can be 
explained naturally within our scheme. First, let us con- 
sider the observed widths of the primary peaks in two 
magnon Raman scattering. We identify two reasons why 
the observed linewidths are very broad, (a) Conven- 
tional spin wave theory yields magnons that do not dis- 
perse along the zone boundary, leading to a singularity in 
the density of magnon states. Consequently, the Raman 
spectrum is dominated by zone boundary magnons. The 
singularity is, however, smoothened by magnon-magnon 
interactions. On the other hand, in the RVB scenario, 
a broad spectrum is obtained even with non-interacting 
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magnons. This is because the low lying mode in Fig. || 
disperses along the zone boundary and there is no sin- 
gularity in the density of states. Thus, the two magnon 
spectrum is not dominated by zone boundary excitations, 
and is broader than the spectrum obtained from linear 
spin wave theory, (b) In addition to this intrinsic width 
of the two magnon spectrum, we expect the line shape to 
be broadened further by the contribution from the spinon 
continuum as is evident from Fig. @. 
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FIG. 3. Density of states of a single spin-flip excitation. 
The dashed line is obtained when only the magnon peak is 
considered. The solid line is 7r -1 S(q,uj). 

To illustrate the above, we calculate the the density of 
states (DOS) corresponding to one spin-flip excitation, 
Y^' S(u) — w(q)). In Fig. the dashed line shows the 
DOS obtained by considering only the magnon peak. To 
incorporate the entire excitation spectrum, i.e., both the 
magnon and the spinon excitations, we also plot (solid 
line) 7T _1 uj). As seen in the figure, the presence 

of spinon excitations leads to a hump at w ~ 2J, fol- 
lowed by a long tail extending to energies u> ~ 4J and 
beyond. These features are clearly absent when only the 
magnon peak is considered. The results indicate that 
the two magnon DOS (relevant to Raman experiments) 
would show a primary peak around lu ~ 3J (from the 
low lying mode) and a broad hump at to ~ 4J (arising 
from the continuum). We expect similar features in the 
mid-infrared absorption spectra. In this case, the spec- 
trum is broader as the two spin-flip excitation can carry 
finite momentum, and the final result involves summing 
over all such momenta To summarize this discus- 

sion, the unconventional dispersion of the magnon in the 
RVB scheme leads to broad primary peaks in the opti- 
cal spectra and the contribution of the spinon continuum 
leads to secondary peaks. 

In conclusion, the observation of dispersing magnon 
peaks in spin 1/2 Mott insulators does not rule out the 
existence of spinons per se. The magnon can be recovered 
as an RVB exciton and the presence of spinons can be 
directly observed by a careful analysis of the line shapes 
obtained from inelastic neutron scattering. We have ob- 
tained, within the framework presented in this paper, the 
characteristic energy scales at which the spinon contin- 
uum can be observed. Dispersing magnons seen in neu- 
tron scattering and the observation of spectral weight 
at high energies in optical experiments, as well as the 



ubiquitous secondary peaks, can all be reconciled by our 
picture. While our results do not prove the existence of 
spinons, they certainly demonstrate how their presence 
modifies the spin excitation spectrum. It is hoped that 
these results would spur further experimental investiga- 
tion of these issues, in particular, a careful analysis of the 
high energy spin excitation spectrum. 
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